Can one extract the electron-phonon-interaction from tunneling data in case of the 

multigap superconductor MgB 2 ? 
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In the present work we calculate the tunneling density of states (DOS) of MgB2 for different 
tunneling directions by directly solving the two-band Eliashberg equations (EE) in the real-axis 
formulation. This procedure reveals the fine structures of the DOS due to the optical phonons. Then 
we show that the numeric inversion of the standard single-band EE (the only available method), 
when applied to the two-band DOS of MgB2, may lead to wrong estimates of the strength of certain 
phonon branches (e.g. the i?2 9 ) in the extracted electron-phonon spectral function a 2 F(u>). The 
fine structures produced by the two-band interaction at energies between 20 and 100 meV turn out 
to be clearly observable only for tunneling along the ab planes, when the extracted a 2 F(uj) contains 
the combination a 2 F^ a (u>)-\-a 2 F^(uj), together with a minor a 2 Fmr(u))-\-(x 2 Fira-(ui) component. 
Only in this case it is possible to extract information on the a-band contribution to the spectral 
functions. For any other tunneling direction, the 7r-band contribution (which does not determine 
the superconducting properties of MgB2) is dominant and almost coincides with the whole a 2 F(u>) 
for tunneling along the c axis. Our results are compared with recent experimental tunneling and 
point-contact data. 

PACS numbers: 74.50.+r;74.20.Fg; 74. 70. Ad Keywords: Eliashberg equations, magnesium diboride. 



There is a growing consensus that superconductiv- 
ity in MgB2 with a critical temperature T c ~ AOK 
(Ref. [I]) is driven by the electron-phonon interac- 
tion (EPI) (for a recent review see Ref. y). An im- 
portant subject to address for a proper understand- 
ing of the surprising physical properties of this mate- 
rial is the character of the order parameter (or super- 
conducting gap): is it constant over the whole Fermi 
surface, or strongly momentum dependent? The idea 
of multiband superconductivity in MgB2******»«*ii is 
supported by many recent experimental results from 
tunneling— point contad*i*ii*«* and specific heat 
capacity measurements'*. These data directly support 
the picture that the superconducting gap has two dif- 
ferent values on two qualitatively different parts of the 
Fermi surface, one A CT for the two quasi-two-dimensional 
a bands and another one A„ for the pair of 3D it bands*— 

While, within first-principles calculations of the elec- 
tronic structure and the EPI in this compound, there is 
an agreementiS on this qualitative picture, still disagree- 
ment is present about the precise values of characteristic 
frequencies and coupling constants. According to most 
calculations*****^ the EPI or, equivalently, the Eliash- 
berg spectral function a 2 F(ui) (EF) is dominated by the 
optical boron bond-stretching i?2 3 phonon branch around 
60 - 70 meV. 

In principle, photoemission or optical measurements 
can also deliver information on the EPl2i, although the 
main experimental tool for the determination of the 
EPI in superconductors so far is the tunneling measure- 
ment. This method has been applied to standard su- 
perconductors with an isotropic (constant in k-space) 



superconducting gap and allowed for the determination 
of the Eliashberg spectral functions in case of many 
conventional low-temperature superconductors (see e.g. 
Ref. HI). 

The spectral EF is obtained from the first derivative 
of the tunneling current 

dlr E 

—— oc Nt — a T Re—F= , (1) 

dv ,/w=*m ^ 

where V is the applied voltage, A(E) is the complex su- 
perconducting gap which depends on energy E, and the 
factor c\t is determined by the properties of the tunneling 
barrier and the corresponding average of the Fermi ve- 
locities of quasiparticles. The standard single-band pro- 
cedure to obtain the EF from the tunneling DOS can 
be found in textbooks* 2 * 2 ^. Another mathematically el- 
egant method has been proposed in Reflzil It has been 
used to investigate conventional (low T c ) as well as high- 
temperature superconductors**. 

Unfortunately, this approach is restricted to momen- 
tum independent s-wave order parameters and cannot be 
used to describe anisotropic superconductors as MgBV 
Nevertheless, there has been a recent attempt to obtain 
the EPI in MgB2 by using this standard approach--. The 
Eig phonon mode has been resolved, but its predomi- 
nance for the electron-phonon coupling was questioned. 
More recently, the E2 g mode has been also resolved in 
point-contact spectra-*. 

The purpose of this paper is to clarify what informa- 
tion can be extracted using this single-band standard 
procedure if one applies it to a two-band superconduc- 
tor. The starting point is the theoretical study of the 
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FIG. 1: The theoretical Eliashberg spectral functions of MgB2 
for the two-band model used in this work (from Ref 6) . The 
thick solid lines correspond to the total a- and 7r-band con- 
tributions. 



quasiparticle tunneling in MgB2-based junctions. The 
superconducting gap functions for the a and 7r band are 
obtained from an extended Eliashberg formalism. The 
parameters for the two-band model utilized in this work, 
which are based on first-principles electronic structure 
calculations^, have been used before for a successful de- 
scription of specific healpi and tunneling'' properties of 
MgB2. The interband and intraband electron-phonon 
spectral functions a 2 Fij(ui), where i,j = it, a (see Fig- 
ure [U and the Coulomb pseudopotential matrix fi*j (see 
Ref l2fJT) are the basic input for the two-band Eliashberg 
theory. The theoretical conductance curves of MgB2 for 
different tunneling directions can be obtained directly by 
solving the corresponding two-band Eliashberg equations 
(EE)1 in the real -axis formulation. The only free param- 
eter is the normalization constant /i in the Coulomb pseu- 
dopotential matrix which is fixed in order to reproduce 
the experimental T c — 39.4 K. 

One may see in Fig. [I] that the ua EPI is dominated 
by the optical boron bond-stretching E2 g phonon mode. 
For other channels there are also important contributions 
from low frequency modes (30-40 meV) and from high 
frequency phonon modes (~ 90 meV). In contrast to the 



FIG. 2: (a) The calculated tunneling DOS in the afc-plane; (b) 
the calculated tunneling DOS along the c-axis. They are both 
obtained by the real-axis solution of the two-band Eliashberg 
equations at T=0 K. The two insets show the fine structures 
of the tunneling DOS due to the electron-phonon interaction. 

case of a conventional junction described by Eq.JTJ, the 
conductance in a MgB2-I-N tunnel junction is a weighted 
sum of the contributions of the DOS of a and ir bands, 
where the weights are determined by the corresponding 
Fermi velocities (plasma frequencies) in the bands and 
by the angle of the tunneling current with respect to the 
ab planed. Figure |21 (a) and (b) shows the calculated 
tunneling conductances in the afr-plane and along the c- 
axis direction. In the case of MgB2 according to Ref0, 

NT(cj) = 0.33JV <r (w) + 0.67A^(w), 

Nj(uj) = O.QUV ff (w) + 0.99N w (u)), (2) 

where A^-(w) and 7V w (w) are the partial superconducting 
DOS. 

The contribution of the it band is always dominant 
even if tunneling is almost in the ab planet In the in- 
sets of Fig. 2 the fine structures due to electron-phonon 
interaction are shown. The maximum amplitude of these 
structures is of the order of 0.5% for measurements along 
the c-axis and 2-3% in the a&-plane. The double-gap fea- 
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FIG. 3: The spectral Eliashberg functions along the afe-plane 
(upper panel) and along the c-axis (lower panel) obtained 
by inversion of the single-band Eliashberg equations (solid 
line) applied to the tunneling DOS of Fig. 2. The dashed lines 
represent the least-square fits of these spectral functions with 
two free parameters for o?F a and a 2 F^ as explained in the 
text. 



tures of the DOS visible in the a&-plane N^ b (uj) should 
be experimentally observable even for a certain amount 
of impurity scattering because the intcrband impurity 
scattering rate appears to be very weakii. On the other 
hand, due to the smallness of the cr-band plasma fre- 
quency in the c direction, the conductance in the c-axis 
direction is almost totally determined by the it band and 
therefore no double-peak structure is expected in the con- 
ductance spectrum. In fact this change in behavior has 
been experimentally observed by spatially well-defined 
tunnel measurements over a step edge in small MgB2 
single crystals 1 ^ and by directional point-contact spec- 
troscopy in larger crystals 17 . 

In order to test the reliability of the numeric inversion 
technique of the standard single-band Eliashberg equa- 
tion when applied to a two-band superconductor, we use 
the calculated tunneling DOS shown in Fig. as an 
input. Of course, in the case of a single-band supercon- 
ductor the inversion should reproduce the shape of the 
single-band Eliashberg function used for the calculation 



of the tunneling DOS. This will not be true anymore 
if one tries to invert a tunneling DOS which is derived 
from a multiband Eliashberg theory. The inverted spec- 
tral functions should correspond to the mixture of the a- 
and 7r-band contributions. We can introduce the follow- 
ing functions: 



l F a (u) = a 2 F arT {u) 
l F n (ui) = a 2 F^{uj) 



a 2 F <y7T (w) 

2 i 



- a 



Namely, these functions determine the normal state 
properties in the two-band model. Tunneling measure- 
ments can only give information on these combinations 
of a Fij(u) — cr, if) which are indicated by thick 
lines in Fig. ^ In order to illustrate this point we show 
in Figure |3 (a) and (b) the results of the inversion of our 
calculated tunneling DOS of Fig. 2 (a) and (b) using a 
standard single-band code (solid lines). The effect of the 
single-band EE inversion is estimated by a least-square 
fit to the inverted spectral functions with the weights 
for a 2 F„ and a 2 F n as free parameters. These fits are 
shown by the dashed curves in Fig. 3. Since tunneling 
along the c-direction -due to the very different plasma 
frequencies in the two bandsi- practically corresponds 
to a single-band case, the inversion properly reproduces 
the a 2 F 7V as one expects. In contrast, in the case of tun- 
neling along the a6-direction (which corresponds to an 
actual multiband situation) the weights are quite differ- 
ent from the theoretical expectation. The results of the 
fit are summarized below: 

a 2 F ab {uj) ~ Ma 2 F a {lu) + 0.16a 2 F„(uj) 
a 2 F c {u) ~ O.Ola 2 F a {uj) + 0.99a 2 F„(uj) (3) 

As one may see the cr-band spectral functions a 2 F acr (uj)-\- 
a 2 F av {uj) play an essential (and amplified) role only if 
the contribution of N a (uj) is significant. The numerical 
simulations show that about 33% of the cr-band contri- 
bution in the tunneling DOS corresponds to a contribu- 
tion of about 66% of a 2 F a (u>) in the effective a 2 F a i,(uj). 
These results are reasonable since A ff = (A W7r + A wcr ) ~ 0.6 
while A,j = (Acrcr + A CT7r ) w 1.23. Somewhat simplifying 
it seems that the contributions of the tunneling DOS to 
the phonon structures are weighted by the correspond- 
ing coupling constants. The inversion of the c-axis case 
results in a a 2 F c (u>) that is almost exactly the sum of 
the a 2 Fi(cu) components taken with the same weights 
present in the sum of the corresponding superconduct- 
ing N v (lj) and N a (u>) (see, Eqs0). In this case, the 
coupling constant from the inversion is almost equal to 
Ajr, and therefore the effective Eliashberg function will 
show strong contributions from low and high frequency 
phonons. 

The above results show that from tunneling measure- 
ments at very low temperatures and by numeric inversion 
of the standard single-band Eliashberg equations one can 
obtain reliable information on the EPI in MgB2 only for 
the 7T band. For doing this we can use the Donetsk's in- 
version program 24 which allows to find the a 2 F(ui) for a 
single-band superconductor. 
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The interesting point is that measuring single crystals 
in the clean limit one should see different phonon contri- 
butions for different tunneling directions, in accordance 
with the two-band model. From tunneling measurements 
exactly along ab and c directions it is possible to ex- 
tract information on a 2 F v (uj) and on a 2 F a (uj), by solv- 
ing the system of equations © and this could be a useful 
method for testing the two-band model and identifying 
the phonon modes responsible for superconductivity in 
MgB 2 . From the experimental point of view this possi- 
bility could remain only virtual since the smallness of the 
phonon structures expected in the tunneling along the c 
axis (see the inset of Fig. 2 (b)) could prevent a correct 
inversion of the EE in the presence of noise. 

In polycrystalline samples the 7r band dominates the 
tunneling current^ and it is only possible to extract in- 
formation about the combination of spectral functions 
a 2 F v (uj), which does not play an important role for the 
superconducting properties of MgB 2 . The recent work of 
D'yachenko et al.^ which claims the experimental de- 
termination of a 2 F[uj), is very likely contaminated by 
strong contributions from the tt band. A value of the 
coupling constant A = 0.9 was reported in Ref^. Un- 
fortunately, such data cannot give information on the 



nature of superconductivity in a two-band superconduc- 
tor as MgB2 which is driven by the interaction in the a 
band. 

The inversion of the calculated afe-plane conductance 
Nj b (uj) shows that a significant contribution to the 
phonon structures comes from a bands and this is re- 
flected in the resulting orF a b(w). Therefore, only in 
junctions with the tunneling current running along the 
ab planes and at low temperature one can observe the 
fine structures of the superconducting DOS produced by 
the electron-phonon interaction at energies between 20 
and 100 meV. In recent point-contact measurements^ 
the anisotropic EPI was observed, though no quanti- 
tative estimate was presented. Further experiments on 
high-quality tunnel junctions are needed in order to ob- 
tain data allowing for a quantitative estimate of the EPI, 
which should also take into account the 'tunneling cone' 
effect, i.e. the distribution of tunneling angles. Ac- 
cording to our result shown in Eq.(J3Jl, separate studies 
of a&-planc and c-axis tunneling conductances may al- 
low a quantitative estimate of the EPI and should thus 
provide a crucial test for the first-principle results of the 
two- band model. 
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